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Abstract 

We argue that the suppression of high pt hadrons discovered recently in heavy ion 
collisions at RHIC may be a consequence of saturation in the Color Glass Condensate. 
We qualitatively and semi-quantitatively describe the data, in particular, the dependence 
upon the number of nucleon participants. We show that if parton saturation sets in at 
sufficiently small energy, then in nucleus-nucleus collisions at RHIC and LHC energies 
the cross sections of semi-hard processes should scale approximately with the number 
of participants, Npart- Our results provide a possible explanation of both the absence 
of apparent jet quenching at SPS energies and its presence at RHIC. Under the same 
assumption we predict that in semi-central and central pA {dA) collisions at collider 
energies the dependence of semi-hard processes on the number of participating nucleons 
of the nucleus will change to ~ (N^rt)^^"^ ■ ^^^^ forthcoming data on dA collisions will 
provide a crucial test of this description. 



The recent results from RHIC 0, |^, ^ suggest that relativistic heavy ion collisions at 
high energies probe QCD in the non-linear regime of high parton density and strong color fields 
B S 0' 0' Pc'-rton dynamics in this new regime is quite different from what is expected in 
perturbative QCD, and is best described by parton saturation ^ resulting in a formation of 
a new type of high density matter, the Colour Glass Condensate (CGC) §, |, 0, 0, 0, 0, |I|]. 
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In this high density regime, the transition amphtudes are dominated not by quantum fluctu- 
ations, but by the configurations of classical field containing large, ~ l/a^, numbers of gluons. 
Even though the coupling as becomes small due to the high density of partons, the fields in- 
teract strongly due to the classical coherence. One thus uncovers new non-linear features of 
QCD, which cannot be investigated in a more traditional perturbative approach. 



This new dynamics appears to be consistent |T^, |I6| with the RHIC data on the mul- 
tiplicity distributions of produced hadrons as a function of centrality, rapidity, and collision 
energy. In the Color Glass Condensate picture, the approximate scaling of hadron multiplic- 
ity with the number of participants Npart, and the logarithmic deviation from it, ~ In Npart 
stemming from the running of the QCD coupling as a function of parton density, has been 



interpreted as a consequence of parton saturation |T5|, |T6[. In this letter we will show that 
QCD saturation gives a very substantial contribution to the suppression of the single particle 
inclusive distributions even at quite large transverse momenta, and can be responsible for (at 
least a part of) the observed approximate Npart scaling of processes at high pf. 



The experiments on high pt hadron production show |T^, |TP| that the inclusive cross 
section is strongly suppressed with respect to the scaling with the number of binary nucleon- 
nucleon collisions Ncoii expected on the basis of the factorization theorem [^] for hard processes 
in perturbative QCD (pQCD). An illuminating way to study this effect is to look at the de- 
pendence of hadron yields in a given pt bin as a function of centrality. The data seem to be 
closer to the approximate scaling with the number of participants, previously established for 
small Pt particles dominating the total multiplicity. It is indeed a surprise that this behavior 
characteristic of soft processes persists up to the transverse momenta (pj ~ 5 -i- 10 GeV), where 
we expect the fragmentation of the incoherently produced jets to dominate. 

The approximate Npart scaling is of course a restatement of the observed earlier suppression 
of the high pt single particle inclusive distributions. One possible explanation of the observed 



suppression is the predicted quenching of the produced jets in hot quark-gluon matter pO|, ^1 
22| . The observation of approximate Npart scaling in the entire range of hadron momenta up to 



~ 6 GeV is however problematic for the scenario in which the jet energy loss is the only reason 
behind the suppression. Indeed, let us assume, following the pQCD approach, that most of the 
hadrons produced at RHIC originate from mini-jet fragmentation. If mini-jet production is a 
hard incoherent process, one expects the initial number of mini-jets to scale with the number of 
binary collisions Ncou- The energy loss in QCD is dominated by the induced gluon radiation and 
therefore as the produced mini-jets lose energy they will emit additional softer gluons. These 
emitted gluons in turn will fragment into hadrons, increasing further the total multiplicity 
of hadrons produced in the event. If the initial number of mini-jets were proportional to 
Ncoii, the measured total multiplicity would then grow even faster with centrality, in a marked 
contradiction to the experimentally observed behavior. This contradiction can be avoided only 
if we assume that the initial production of moderate pt mini-jets, the fragmentation of which 
gives the dominant contribution to the measured multiplicity, is suppressed with respect to 
Ncoii scaling. (The jets with really high pt because of the small production cross section do 
not contribute to the total hadron multiplicity in any significant way, and their production 
according to Ncoii scaling, and subsequent energy loss, would not lead to any confiict with the 
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observed scaling of the total multiplicity). 

The main goal of our letter is to show that the approximate Npart scaling of mini-jet produc- 
tion at moderately high pt (up to ~ 6 ^ 8 GeV at RHIC energies) has a natural explanation 
in terms of parton saturation and the Color Glass Condensate. We also demonstrate that the 
Npart scaling exhibits very non-trivial properties of dense partonic systems not only at momen- 
tum transfers around the saturation scale Qs-, but also at much higher momenta, naively as 
large as ~ Q^/Aqcd- For heavy nuclei, this new momentum scale is 0(5 — 10 GeV) at RHIC 
energies, and even larger at the energies of the LHC, (9(25 — 50 GeV). We will in fact find 
that the correct scale for nuclei is somewhat reduced relative to this estimate, but is still large 
compared to Qsat- We find this scale to be about 4 GeV at RHIC and of order 10 GeV at the 
LHC. 

Our arguments are based on the following foundations. First, in the Color Glass Condensate, 
the saturated parton density provides a new dimensionful scale; as a consequence, the structure 
functions at momentum transfer depend on a single variable r = Ql/Q^ ("geometrical 
scaling"). This has been proven ^ ^ ^ ^ ^ ^ ^ for Q < Qs- Moreover, it has been 



observed that the deep inelastic scattering data at small x < 0.01 exhibit geometrical scaling 
in the entire kinematical region of accessible at HERA, and it has been shown that 
the scaling indeed holds in a wide range of Ql < < Q^/Ag^^. Second, the properties of 
QCD evolution are modified in the domain of high parton density, where geometrical scaling 
holds. As we will discuss below, the anomalous dimension of the gluon density in this region 
is approximately equal to 1/2, as was computed in [|, |2^, Third, in the range of 



< < Qs/Aqq^ the saturation scale for nuclear targets scales according to Ql ~ A^^^, a 
behavior which was argued to be responsible for the centrality dependence of hadron multiplicity 
at RHIC 0, |T§. 



Let us begin by considering an external probe with virtuality interacting with a target 
of longitudinal size L and an average (3-dimensional) parton density n. The probability of 
interaction in the target is determined by the ratio n = L/X of the target size to the mean 
free path of the probe A = (crn)~^, where a ~ Ois{Q'^)/Q'^ is the cross section of the probe 
scattering off a parton. At small Bjorken x, when the coherence length ~ l/{'mx) (m is the 
nucleon mass) exceeds the target size, the probe can interact coherently with partons located 
at different longitudinal coordinates, and the only relevant parameter characterizing the target 
becomes the density of partons in the transverse plane, p = nL; for a nuclear target, this implies 
p ~ A^^^. The interaction probability is then determined by a dimensionless "parton packing 
factor" K ~ p{x) cxsiQ"^) / Q'^ ■ The boundary of the saturation region, and the corresponding 
value of "saturation scale" Ql{x), are given by the line in the (x, Q^) plane along which the 
interaction probability is of order one, k{x, Q"^) ^ 1. 

We will now write the packing factor in the form WM suited for nucleus-nucleus collisions. 



where a centrahty cut can be used to select events with different transverse densities ppart of 
nucleons participating in the collision ("participants"): 

. = 'j0^^GM-,^ (1) 
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where xG{x, Ql) is the gluon structure function; we assume that x is sufficiently small for 
gluons to dominate the parton densities. Consider the kinematical regime of small x and high 
in which both as InQ^ and as ln(l/x) are large; in this domain both logarithms have to be 
resummed, and the solution to DGLAP evolution equations for xG can be written in the 



"double log approximation" (DLA), which allows us to re- write Eq. (|Ij) in the following naive 
form, neglecting the influence of saturation effects on the QCD evolution: 



exp 



part 



(2) 



we have assumed that at the initial point Xq of the evolution in x the transverse momentum 
squared of partons is Qg, and have introduced the following notations: y is rapidity, with 
y - yo = ln(xo/x), ^ = ln{Q'^/Ql), ^part = ln{ppart/ Po) ~ |lny4, A is the atomic number, 
Po is parton density in a nucleon, and ds = {.Nc/'n)as in the large approximation. 

We now investigate the effect of saturation on the function k. Our treatment so far assumed, 
in the spirit of pQCD, that evolution starts at some initial fixed scale Qq ~ ^qcd separating 
perturbative and non-perturbative domains, and that ^ Qg- Traditionally, one assumes 
that the dynamics at the scale and below is driven by non-perturbative phenomena and 
should be described by a priori unspecified initial conditions for the structure functions, which 
are then constrained by experimental data. The parton saturation approach is different in that 
it provides theoretical tools to compute the parton distributions once their density becomes large 
enough. If saturation occurs in the nuclear wave function at some Xq and at the corresponding 
scale Q1{xq), in performing QCD evolution of parton densities to smaller values of x and larger 
values of we should therefore provide the saturated Color Glass Condensate distribution as 
the initial condition. 

It is convenient to recall the well-known property of DGLAP equation |^ (see Ref. |Q for 
example) which allows to write its solutions in a form expressing their Green's function-like 
properties: 

= f de niy -yo,^- e) <yo, O (3) 

Using the DLA form of the solution we see that the integral in (|]) has a saddle point at 
^' = ihjo/y)- If we consider y ^ yo {x xq), then the dominant contribution to the integral 
comes thus from ^' — »• 0, corresponding to the initial condition at Qq <^ Q^. If, however, 
Xq is sufficiently small for the saturation to set in, the initial condition will be provided by 
McLerran-Venugopalan model at the scale C,' = C,sat{yo) = ^^{Qliyo) / Qq) ■ The integral over 
^' will then be dominated by ^' ~ C,sat{yo)- This results in a different solution, which instead 
of Eq. (0) is now given by 



exp 



4aS (y - yo) - ^satiyo)) - ^ + ^satiVo) 



(4) 



The saturation boundary can be found from the condition K = 1. This leads to the following 
equation for the saturation scale at rapidity y: 

4 asi^satiy) - ^satiyo) ) " ^/o) = {^sat{y) " ^satiyo)f ■ (5) 
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This equation has a simple solution 



^sat{y) = 4:as {y - yo) + 6at(l/o) • (6) 

Since Q^xq) is proportional to at relatively low energies, when y ^ yo (see Eq. (||) or the 
McLerran - Venugopalan formula [§)), one can still conclude that Ql{y) oc As. For nuclear 
collisions, this means that the product of the nuclear overlap area Sa and the saturation scale 
is proportional to the number of participants, SaQI oc Npart HI, HH- 

Let us establish the leading behavior of the solution at high . At very high ^ 3> 
'^(.satiy), the solution tends to 



Q 



^ - ^> (7) 



Q 



modified by the double logarithmic factor of exp (^^J4as{y — ?/o)'C)- This result does not exhibit 
geometrical scahng, and is the same as the naive perturbative result. Consider, however, the 
region where 

^-Uy)^Uy)-Uyo), (8) 

which corresponds to a vast kinematical domain of 

^ Q'siyoY ^ ^ 

note that this is precisely the condition for the observed geometrical scaling to hold in deep- 
inelastic scattering ||2^. In this domain, we should expand the exponent of Eq.(0) in powers of 
the ratio AC, / {C,sat{y) —C,sat{yo)) ^ 1, where we have introduced the notation = C,—Csat{y) = 
ln[Q'^/Ql). The expansion to first order yields 



4^5 - Z/o) - ^satiyo)) - ^ + isatiy^ 



= {Ui{y) - Ut{y^)) y 1 + ^ ^^y^^}^ - - iUt{y) - Ut{yo)) ^ -\ A^, (lo) 

where we have used Eq.(|I). This result translates into the following dependence of the 
packing factor: 

.^(^^)^ 

The appearance of the square root clearly shows that the properties of QCD evolution are 
dramatically affected by the saturation phenomena in the entire region (^, much more broad 
than the originally anticipated ~ Q^{y)- The comparison of and (|TT|) reveals that in the 
presence of Color Glass Condensate the gluon density acquires the anomalous dimension equal 
to 1/2. 

Let us now turn to the discussion of semi-hard processes in hadron and nuclear collisions. 
Due to the AGK cancellation [^], the inclusive cross section for large pt jet production is 
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Figure 1: The simplest Mueller 
diagram for the high pt jet in- 
clusive production. 



described by the Mueller diagram shown in Fig. It can be written as (see Refs. |3J, qE 



^^p = wh^, I dkl as^M{y -y.kl)vM{yAp-k)l\ (12) 



where (/j's are the unintegrated parton distributions related to the packing factor k discussed 
above by simple relation ip ~ Sai^a/(^s^ and Y is the rapidity of the measured jet. Using the 
formula Eq.(|lT]) for ip ~ SaK'a/^^s in the domain (|]), and noting that in this case = p^, we 
see that Eq. (|12D leads to 

oc SaQI/pI - N,^^,/pI (13) 



dyd?pt Sa d?P 

where we have used that Npart oc QISa- Note that the number of produced hard particles 
scales with the number of participants. The Npart scaling given by Eq. (|T3|) can be formulated 
in the following way: 

1 dN . , . , , „ „ Qliy) ^ / N 

IS almost mdependent of centrality cut tor — - — > Pt > Qs[y) , 



Npart dyd^pt A 

where A = Qs{yo) is the saturation scale at some yo < y marking the onset of saturation. This 
is the central result of our paper. 

Of course, "genuine" hard processes still scale with the number of collisions. Indeed, consider 
Pt > Qll ^] then according to Eq. (0), Ql/pj and the formula Eq. (^) yields the number 

of hard particles proportional to SaQ^pI ~ ^coii/Pt iii accord with the pQCD factorization 
theorem. 

We have demonstrated that at high energies corresponding to y ^ yo the effects of the 
Color Glass Condensate persist up to very large values of pt; how far do these effects extend at 
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RHIC and LHC energies? The answer depends on the energy at which the saturation effects 
set in. It seems reasonable to assume, on the basis of the analysis of centrality and energy 
dependence of hadron multiplicities [iZ, 15, 16], that the minimal value of saturation scale is 



about Qliyo) ^ 0.5 GeV^. In central Au—Au coUisions at RHIC energy, where Ql{y) ^ 2 GeV^, 
the boundary of the kinematical domain fully dominated by saturation extends thus up to 
~ V2 Q^siy)/Qs{yo) - 4 GeV. At LHC, where ^ 4.5 GeV^ [|16l, this boundary extends 



up to pt — 9 GeV. It is also important to note that in peripheral collisions, as Ql{y) becomes 
small and approaches the value of Q'^{yo), we return to the familiar situation in which the 
boundary of hard processes sets in at a rather small value of ~ 1 GeV, and the number of 
produced mini-jets, and hadrons, thus scales proportionally to the number of binary collisions. 
However, as the centrality increases, and Ql{y) grows, the number of the produced semi-hard 
particles begins to deviate from the binary collision scaling; the transition to the Npart scaling 
occurs. 

While these arguments are well defined and grounded theoretically, reliable computations 
will require further improvement of the quantitative understanding of saturated parton distri- 
butions. Nevertheless, we will present a calculation of the semi-hard particle production using 
a simple input distribution stemming from Eq. @): 



¥^A{y,Pt) = exp 



^4as {y-yo)m)] ^A{yo,Pty, (14) 



where ^{pf) = ln(Pt /Q^(l/o))- To bring this formula in line with the standard description 
of the deep-inelastic scattering (DIS) data, we have to modify it by slightly changing the 
anomalous dimension. Indeed, it has been demonstrated |^ that the following simple form of 
the anomalous dimension in the Mellin moment variable u (in the leading order in a^): 

7M = «5 (^-l) (15) 

describes the DIS data well. This "phenomenological" anomalous dimension leads to a multi- 
plicative factor of (Qs/Q^)""^ formula in the x-pt representation. The resulting anomalous 
dimension is consistent with the result found from analysis of the BFKL equation and satu- 



ration in leading and next to leading order |40, O, E2l H4[. We use the running coupling 



in (p!5|), freezing its value at small virtualities at a™"^ = 0.45; therefore, at very high ener- 
gies, when as < d^°-^ in the entire region (^), the additional multiplicative factor turns into a 
constant, and we return to the anomalous dimension of 1/2. 

For the function ^\{yQ-,pX)-, describing parton distributions in the Color Glass Condensate 
at the onset of saturation, we will use the formula 



as Ql{x) 

<iM=Pl±i),„i^S±i - 1. (16) 

^ ' V4r + 1 V4r + 1 - 1 ^ ^ 
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(1/N^^^) d^N/d^p^dr|(GeV-^) 
Au + Au W=130 GeV at r[= 




Figure 2: The pj-spectrum of charged hadrons in central (0 — 10% centrahty cut) Au — Au 
coUisions at -v/i = 130 GeV. The data is from lll^ . 



where Sa is the nuclear overlap area. At large pt, this distribution falls off as 1/p^, corresponding 
to the classical bremsstrahlung, and at small "pt it has the behavior ~ ln(Qg/p^), established 
previously for the Color Glass Condensate |^6|, It agrees well with the numerical solution 
of non-linear evolution equations 



26| 



its properties will be presented elsewhere ||32| 



the derivation of Eq. ( [T6| ) and detailed discussion of 
The value of the initial rapidity should be 
chosen to correspond to the onset of coherent partonic phenomena; we choose = yi, — ln{l/xo) 
{yb is the beam momentum) corresponding to the value of Bjorken xq ~ 0.1; this marks the 
beginning of the regime in which the coherence length ~ l/(mxo) starts to exceed the typical 
inter-nucleon distance. 

The cross section of hadron production can then be obtained by convoluting the result of 
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Figure 3: Centrality dependence of hadron yields per participant pair in Au — Au collisions at 
^/s = 130 and 200 GeV in different pt bins; the yields are normalized to the yield in peripheral 
collisions. Upper solid lines show the behavior expected in perturbative QCD. 
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Figure 4: Same as in Fig.|^ in the pt bins of 6 and 9 GeV. 
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calculation according to Eqs.(0,0) with the jet fragmentation function Dfrag{z,pt)'- 

= J "^'d^? " ; (17) 

we also introduce in the distribution functions the common ~ (1 — x)^ factors enforcing the 
behavior at x — 1 prescribed by the quark counting rules. 

The result for the charged hadron spectrum in central Au — Au collisions at =130 
GeV is compared in Fig.|^ to the experimental data from PHENIX Collaboration. While the 
agreement is not perfect, the calculated cross section does agree with the data within error bars 
(only statistical error bars are shown), which is remarkable in spite of the crude approximations 
we made and our use of the standard fragmentation functions fitted for the use with pQCD 
calculations. In Fig. |^, we show also the dependence of particle yields in different pt bins on 
centrality. While the result admittedly depends on the approximations we made and is not very 
robust, it does show that up to quite large values of pt the yield of particles deviates strongly 
from the scaling with binary collisions, shown by the upper curve on Fig.|^. 

Further work is clearly needed before a reliable calculation of high pt cross sections in the 
Color Glass Condensate approach can be performed. Nevertheless, it is already clear that high 
density QCD effects have to be taken into account in the interpretation of the discovered at 
RHIC suppression of high pt particles. 



The mechanism proposed here may also explain the old puzzle the absence of apparent 
jet quenching in the SPS data on high pt neutral pion production. Indeed, at relatively 



small energy of the SPS (a/s ~ 17 GeV) the transverse momentum spectra of the produced 
particles are rather steep, and even a small jet energy loss would induce large deviations from 
the perturbative scaling with the number of collisions. These deviations were not observed: the 



data appear consistent with the perturbative QCD calculations At the same time, many 
of the "soft" properties of the collisions, as well as the suppression of charmonium, point to the 
presence of strong collective effects (see, e.g., |^ ). 



Basing on the fact that the prediction of KLN model |jT6[ for centrality dependence of hadron 
multiplicities in the ^/s = 20 GeV run at RHIC appeared to be successful , it may be possible 
to assume that saturation sets in around SPS energyQ. Since saturation provides favorable initial 
conditions for thermalization the emergence of collective effects in "soft" observables 

and charmonium suppression would thus be likely. However, since for y ^ yo the domain (|^) 
Qliy) < < Qtiy) / Qliyo) shrinks to zero, the hard processes above the saturation scale 
exhibit the usual perturbative behavior, and no apparent jet quenching should appear. 

The effects of our suppression mechanism begin to fail around 6 GeV for the 130 GeV 
data and somewhat higher for the 200 GeV data. We might have been able to work harder 
and find some marginally acceptable parameterization of the gluon structure functions which 
would extend the agreement to larger pt, but this would not be very natural. Part of the 
problem is at these large pt values, we become sensitive to the gluon structure function within 

*We do not believe that extrapolation to still lower energies would make sense since the coherence length 
becomes much smaller than nuclear size in that case. 
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the nuclear fragmentation region x > 0.1, and this is a region which is more comphcated than 
that for which our naive saturation model can accommodate. It is also true that there must 
be some effects of jet energy loss in the media produced in heavy ion collisions, and the degree 
to which our saturation mechanism is dominant or subdominant relative to this mechanism is 
impossible for us to assess a priori. 

The dA data at RHIC will allow to disentangle the effects of jet energy loss from gluon 
saturation. Saturation is a property of the nuclear wavefunction, and in so far as it modifies 
jet production in AA collisions, there will be a corresponding effect in dA collisions. Let us 
discuss this in more detail. In dA collisions at RHIC energy, with the exception of peripheral 
interactions where the density of participants, and thus the saturation scale, is too small, the 
gluon distribution in the nucleus A is saturated and in the domain Qg < Q < Q^/ A is thus 
given by 

SaI^A Sa ( QXa\ /-ion 

The deuteron distribution is given by ipd ~ Sdi^d/oig] our conclusion on the scaling with the 
number of participants from the nucleus A will be independent on the explicit assumption 
about ifd- The dependence of semi-^hard processes in the saturation region on the number of 
participants from the nucleus A N^^^^ will now be given by 

^ = Ta^I^ ~ ^^-^^ ~ i^--) • ^^^^ 



where we have used that in dA collisions a ^ ^part- We predict the behavior ([l9| ) for semi- 



hard processes to hold in dA collisions at RHIC with the exception of peripheral collisions, where 
the saturation scale ^ is too small. The yield of particles per participant from the nucleus A 
in these collisions is thus predicted to decrease as {Np^^^)~^/'^, deviating from the scaling with 
the number of collisions expected in perturbative QCD. We expect the saturation effects to set 
in around N^^^ ~ 6, corresponding to the impact parameter of dAu collision 6 = 5 -i- 6 fm 



We thus predict that around N^^^ ~ 6 the yields of high pt particles would begin to deviate 



from the scaling with the number of collisions Ncoii ~ N^^^; the yield per participant will start 



part 

■ collisions Ncoll ~ Apart) 

to decrease as (Aj^"J~^/^. In 15% most central dAu events, where N^^^ ~ 12 |5^, we therefore 
expect to see the normalized yield of (6/12)^/^ ~ 0.7, corresponding to ~ 30% suppression 
of high Pt particlesQ. The scaling of semi-hard processes in dAu collisions with centrality at 
RHIC energy can thus expected to be drastically different from that observed previously at 
fixed target energies, where at high pt the yields of particles were proportional to the number 
of collisions. 

On the other hand, jet energy loss is expected to be a minor effect for dA collisions since the 
jets are produced in a region far from the fragments of the nucleus for recent predictions, 
see |5l. 



^Numerical calculations according to the formulae given above give somewhat smaller, but close, ~ 25% 
suppression effect. 
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In addition, the dA collisions provide the opportunity of studying jet production in the 
fragmentation region of the deuteron, and this probes much smaller values of x of the nucleus 
than is the case for centrally produced jets. One expects the saturation scale to grow as 
EH Ql A^y) = Qs a(0) 6xp(A?/), with A ~ 0.25. At 3 units of rapidity away from the central 
rapidity region toward the deuteron fragmentation region in dA collisions at ^/s = 200 GeV, 
the saturation scale thus grows to about Ql = 3) ~ 4.5 GeV^. The approximations we 
are forced to make for centrally produced jets in AA collisions should be under much better 
control for jets arising in the fragmentation region of the deuteron; recently, this problem was 



addressed in 56 



n, mi, [n, IT 
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